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Abstract. Dynamics of stacked intrinsic Josephson junctions (IJJ) in the high-Tc 
superconductors is theoretically investigated. We calculate the current-voltage characteristics 
(CVC) of IJJ and study the breakpoint region on the outermost branch of the CVC for the 
stacks with 9 IJJ. A method for investigation of the fine structure in CVC of IJJ based on 
the recording the "phase-charge" diagrams is suggested. It is demonstrated that this method 
reflects the main features of the breakpoint region. 



1. Introduction and model 

Study of intrinsic Josephson junctions in HTSc, like I^SV^CaC^Os-HS shows interesting 
physical features. In Refs.[IJ [2 [3] we studied the multiple branch structure of the CVC of 
IJJ and showed that the branches have a breakpoint (BP) and some breakpoint region (BPR) 
before transition to the another branch. The BP is determined by the creation of the longitudinal 
plasma wave (LPW) with a definite wave number k, which depends on the coupling parameter 
a, dissipation parameter /3, the number of junctions in the stack, and the boundary conditions. 

Here we show that using the phase-charge diagrams (phase portraits) give us an additional 
method for investigation of the fine structure in the CVC. We use the (CCJJ + DC) - model 
[I] to investigate the phase dynamics and CVC of coupled system of Josephson junctions. The 
system of equations in this model has a form 

-jjpfl = [I ~ sin ipi - P—j^) + «(sm <pi+i + sin <^_i 
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Here ipi is the gauge-invariant phase differences ipi(t) = 9i + i(t)—0i(t) — jr ff +1 dzA z (z, t) between 
superconducting layers (5-layers), 6i is the phase of the order parameter in S-layer I, A z is the 
vector potential in the barrier. The current, voltage and time are normalized to the critical 
current I c , Vq = and inverse of plasma frequency cj" 1 , respectively. The details of simulation 
are available at Refs. [5j[3l[6]. We calculate the charge on S-layer as Qi(t) = a(Vi(t) — V/_i(t)), 
where V[(t) is the voltage drop on insulator layer 1 at time moment t. We start from the 
current value 1=1.2 and decrease it up to the transition to another branch. We record the time 



dependence of the phase difference between layers I and I — 1 and the electric charge on the layer 
I at some selective points of the BPR. 
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Figure 1. (Color online) The BPR part of the outermost branch for the stack with 9 IJJ at 
a = 1, f3 = 0.2 and periodic boundary conditions. The symbol V indicates the points at which 
the time dependence of the phase difference and the charge on the layer are investigated. 



2. Phase-charge diagrams 

The BPR part of the outermost branch of CVC for the stack with 9 IJJ at a = 1 and (3 = 0.2 
and periodic boundary conditions is presented in Fig. [TJ It indicates the current values, at which 
we have investigated the phase-charge diagrams. The phase-charge diagram is a trajectory of 
junction in its phase space: (ip, Q). We suggest to use this diagrams as a tool to probe the 
dynamics of system. For periodic systems like simple harmonic oscillator, the phase portraits in 
(<p,<p) space are circles with their centers at the origin of coordinate system. For the system of 
IJJ it is natural to investigate the phase-charge diagrams. 

Fig. [2] shows such phase-charge diagrams: the variation in time of the modular of charge on 
the first superconducting layer Qi and phase difference ip\ in the first Josephson junction of the 
stack. In these figures the radial axis shows the charge Q\ and azimuthal is the phase ip\. 

The beginning part of BPR corresponds to the onset of the parametric resonance: charge on 
layers is growing. In this region, particularly, / = 0.57126 (Fig. [2k.) the Josephson frequency is 
still high than LPW frequency, so the trajectory in phase-charge space has some small thickness. 
At I = 0.569 (see Fig. [2b) the conditions u j = 2ulpw is practically ideally fulfilled, charge on 
layers has reached its maximum amplitude, and thickness of trajectory in phase-charge space is 
just thickness of line. Results of FFT analysis (not presented here) show that at this value of 
current we observe the peaks corresponding to the LPW only and Josephson frequency does not 
manifests itself. The system follows this trajectory all time. As we show below, the results of 
analysis of autocorrelation function, presented in Fig. [3] confirm this interpretation [6]. 




Figure 2. (Color online) The phase-charge diagram (the charge on the first superconducting 
layer versus the phase difference in the first junction) at different value of bias current. The 
radial axis is a charge, the azimuthal axis is a phase difference. 



At I = 0.566 (Fig. [2}:) the trajectory in phase-charge space is wide again. The resonance is 
passed: ujj > 2ujlpw- Here the amplitude of charge on the layers is merely constant but there 
is a beating of LPW in the stack. A new large period of system appears [7j, but the system 
after some cycles repeats the same trajectory. Analysis of charge correlations shows that the 
behavior of system is regular, while Fig. [2ti corresponds to a chaotic behavior of the system. 

There is a correspondence between phase-charge diagrams and autocorrelation functions, 
presented in FigJ3J The smallest value of r = tq with C a (r = tq) = C a (r = 0) is a period 
of system. If the autocorrelation function has not such points, then the trajectory of system 
in phase-charge diagram is open. There is such r = tq in the cases presented in Figs. [3b and 
[3b which satisfies the mentioned condition and that's why in Figs [2b and [2b the trajectories 
in phase space are closed. But in Fig. [3U the autocorrelation function doesn't show such tq 
and hence Fig. [2H show open trajectories. Behavior of the system in this case is chaotic. As 
we mentioned above the thickness of trajectory in Fig[3b is related to the growing charge on 
the superconducting layers and this is a reason of the amplitude's decreasing of autocorrelation 
function with time. Note that to show the very small decreasing of autocorrelation we zoomed 
the region in vertical axis of FigJSb,. The value r is introduced in the definition of autocorrelation 
function 

Cf (r) = lim — 1 Qt(t - r)Qi{t)dt (2) 



C a (T = 0) is a maximum and its value is C?(t = 0) =< Qf(t) >. At I = 0.566 (Fig. Eb) the 
trajectory in phase-charge space is wide again. The resonance is passed: ojj > 2ui lpw- Here 
the amplitude of charge on the layers is merely constant. 
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Figure 3. (Color online) The charge autocorrelation function C a for the first layer at different 
values of bias current in the breakpoint region. 

3. Summary 

We found that the phase-charge diagram reflects the dynamical behavior of IJJ. It was shown 
that the results of the phase-charge diagram analysis are in agreement with the results of the 
autocorrelation function analysis. We demonstrated that using the phase-charge diagrams (phase 
portraits) give us a powerful method for investigation of the fine structure in the CVC. 
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